
VIBRAÇAO TRANSIENTE 





4.1 INTRODUÇÃO 

Quando um sistema dinâmico c excitado pela aplicação súbita de uma excitação 
F{t) não-periódica, tal como a representada na Fig. 4.1-1, a resposta a este tipo 
de excitação c denominada resposta transiente, uma vez que não são gcralmentc 
produzidas oscilações de estado permanente. Tais oscilações ocorrem na frcqiiôncia 
natural do sistema, variando a amplitude de uma maneira dependente do tipo da 
excitação. 

Inicialmcnte, estudamos a resposta de um sistema mola-massa a uma excitação 
de impulso, por ser este caso importante para a compreensão do mais geral problema 
de transientes. 


4.2 EXCITAÇÃO DE IMPULSO 

Impulso é o tempo integral da força, c o designamos pela notação /■’ 

l- = \TU)<lt ■ (4.2-1) 

Encontramos comumente uma força, de muito grande magnitude, que atua durante 
um período de tempo muito curto, mas com um tempo integral que c finito. Essas 
forças são denominadas impulsivas. 



Figura 4.1-1 Excitação não-pcriódica. 


A Fig. 4.2-1 mostra uma força impulsiva de magnitude Fje com uma duração 
de e. À medida que e se aproxima de zero, tais forças tendem para infinito. Entre¬ 
tanto, o impulso definido por seu tempo integral c F, que é considerado finito. 





Figura 4.2-1. 

Quando fi é igual à unidade, tal força no caso restrito de c -> 0 c denominada 
unidade de impulso ou a função delta. A função delta quando r = f é identi¬ 
ficada pelo símbolo 5(í-|) e tem as soguinte.s propriedades 

5(f - I) =‘0 para todos os valores de t ^ 

= maior que qualquer valor admitido para t = ^ (4.2-2) 

r<5(r - Oí/r = 1,0 0<í'<üo 

J 0 

Se 5 (/ — Q é multiplicada por qualquer função de tempo f{t), como indicado 
na Fig. 4.2-2, o produto será sempre zero, exceto quando í = ?, c sua integral 
será 


Desde que Fdt = mdv, o impulso F atuando sobre a massa resultará numa sú¬ 
bita mudança na sua velocidade igual a F/m sem apreciável mudança no seu desloca¬ 
mento. Quando da vibração livre, constatamos que o sistema mola-massa não-amor- 
tecido com condições iniciais x(0) e x(0) comporfava-se de acordo cora a equação 

... á:( 0) c<»n t -.L -yfCW nr\c /'A f 


X rrc -i^sen coj -h x(0) cos coj 
CO,. 


Por isso, a resposta de um sistema mola-massa inicialmcnte em repouso e excitado 
por um impulso F é 


Quando o amortecimento está presente, podemos iniciar com a equação de 
vibração livre 

.V ^ - ,V't’‘'“"'sen(.yi — cuj — (f>) 


.V ^ - -V't’‘'“"'sen(.yi — cej — (p) 
c substituindo as condições iniciais acima, chegamos à equação 

.V - -sen, (4.2-5) 

/nor„ y I ■ - 

A resposta ao impulsp unitário c importante para os problemas de transientes, 
e c identificada pela designação especial siO- Nestas condições, quer se trate de um 
caso amortecido óu não-amortecido, a equação para a resposta impulsiva pode ser 
expressa como se segue 


.V - Fg{t) 


onde o lado direito da equação c dado por uma ou outra das Eqs. (4.2-4) ou (4.2-5). 


4.3 EXCITAÇÃO ARBITRÁRIA 

Tendo a resposta g;(í) para um impulso unitário 'de excitação, é possível estabelecer 
a equação para a resposta do sistema excitado por uma força arbitrária f{t). Para 
este desenvolvimento, consideramos a força arbitrária como sendo uma série de im¬ 
pulsos, conforme a Fig. 4.3-1. Se examinamos um dos impulsos (o que está hachu- 
rado) no tempo t — Ç sua força é 
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/' = /(í) Aí 

e sua contribuição para a resposta no tempo / depende do tempo decorrido (t - J), 
ou 

/(í)Aíg(t-í) 

Sendo linear o sistema que estamos considerando, o princípio dc superposição pre¬ 
valece. Desta forma, combinando todas essas contribuições, a resposta para 



Figura 4.3-1. 


. a excitação arbitrária f{t) é representada pela integral 

A-(/) - r/(í)í(í ■ «í)'/í (4.3-1) 

j 0 

A integral acima é chamada dc Convotução integral ou algumas vezes referida como 
a superposição integral. 

Estabelecendo r = (f - ?) encontramos outra forma desta equação. En¬ 

tão, i = t - T, d% = -dr, e obtemos 

a(í) -- í f(l - ■t)g{x)dx (4.3-2) 

0 

Quando t é maior que a duraçãp do pulso, isto c, o limite superior da 
Eq. (4.3-1) permanece cm t^, porque, então, a integral pode ser expressa corno 


A(í) - J'' /(í)g(/ - í) r/í -h /(Off(í - ■ í) r/í 

- í"/(Í)g(í - í) r/í, í>í, 

j 0 

Aqui, a segunda integral c zero, uma vez que /(?) = 0 para ? > íj 
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(4.3-3) 


Excitação base. Muitas vezes o suporte do sistema dinâmico c sujeito a um movi¬ 
mento repentino estabelecido por seu deslocamento, velocidade, ou aceleração. A 
equação dc movimento pode então ser expressa em termos do deslocamento relativo 
z = X - y como se segue 

ê 2íco„i -1- co^z = ~y (4.3-4) 

e daí, todos os resultados para o sistema excitado-força aplicam-se para z no sistema 
excitado-base, quando o termo A’o/m é substituído por -jt ou o negativo da ace¬ 
leração de base. 

No caso de um sistema não-amortecido, inicialmcnte em repouso, a solução 
para o deslocamento relativo torna-se 

‘ (4.3-5) 

n J 0 

Exemplo 4.3-1 

Determinar a resposta dc um sistema de um grau dc liberdade à excitação 
degrau representada na Fig. 4.3-2. 



Figura 4.3-2. Excitação da função degrau. 


Solução: Considerando o sistema nao-amortccido, temos 

gl/) — —sen ctí,/ 
inO}„ 

Substituindo na Eq. (4.3-1) a resposta do sistema não-amortecido é 
.v(/) : j scnco„(/ ■ íjí/í 

'"W„ j „ " ^ 

- ^( I - cos mj) 

Este resultado indica que a resposta máxima à excitação degrau de magnitude Fq é 
igual a duas vezes a deflexão estática. 


Para um sistema amortecido, o proceçso é repetido com 

,?(t) =--- ■ sen Vl C* (o„i 

mcu„,v'l -í- 

ou, alternativamente, podemos considerar çimplesmcntc a equação diferenciai 
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X + 2^a>„x -h colx = ^ 
m 

cuja solução é a soma das soluções da equação homogênea e da solução particular, 
a qual para este caso é . Assim, a equação 


m(ül 


ajustada às condições iniciais de Ji:(0) = :í(0) = 0 resultará na solução que é dada 
como 



onde 



A Fig. 4.3-3 mostra um gráfico de xkjFa versus u>„t, com f como parâmetro, e é 
evidente que a resposta máxima é menor que 2Fo/k, quando o amortecimento 
está presente. 



Exemplo 4.3-2 

Considerar um sistema mola-massa não-amortecido onde o movimento da 
base é especificado por uma velocidade de pulso da forma 

>■'(0 " 

que é representada na Fig. 4.3-4 juntamente com sua taxa de variação de tem¬ 
po a = f . 

Solução: A velocidade de pulso para f = 0 dá um salto repentino de zero pára 
Vg e sua taxa de mudança (ou aceleração) é infinita. Admitindo que /a t/t = v, 
a súbita mudança na velocidade para r =. 0 é satisfeita por 

f õ(i)(l! r v„ 

0 

Assim, a aceleração da base torna-se 
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Figura 4.3-4. 


Da substituição de y{t) na Eq. (4.3-5) resulta 

<t) j' (^(í) - -ie-' '•[ sen,«„(/ - í) t/í 

-- |—— (4>Jo sencu,/ — cos coj] 

Exemplo 4.3-3 

Uma massa m fixada a uma mola de rigidez k é sujeita ao impulso repetido 
de duração desprezível, a intervalos de t,-, conforme a Fig. 4.3-5. Deter¬ 
minar a resposta de estado permanente. 

Solução: Entre cada dois impulsos, o sistema está em vibração livre na sua frcqüên- 
cia natural = -s/k/m. Fazendo t -■ 0 imediatamente após cada impulso. 



Figura 4.3-5. Impulso repelido num sistfima mola-massa. 


O deslocamento c a velocidade podem ser expressos como 

X = A sen (caj -|- (a) 

X — a>„A COS (cúj -f 0) (b) 





e em conseqücncia, com t = 0, temos 

x(0) ^ A sen 0 (c) 

Á(0) ~ 0}^A COS 0 (d) 

Justamente antes do próximo impulso, o deslocamento c a velocidade são 

.v(t,) = A sen -j- 0) (e) 

í(t,.) '= a„A COS (a)„Ti 1- 0) (0 

onde T,- é o intervalo de tempo entre os impulsos. O impulso atuando nc.sfc tempo 
aumenta a velocidade subitamente - de Fim, embora o deslocamento permaneça 
essencialmente sem modificação. 

Se é atingido o estado permanente, repetem-se o deslocamento e a velocidade 
depois de cada ciclo. Deste modo, podemos escrever 

/I sen 0 = /( sen {<o„x, 0) (g) 

(ú„A COS 0 = (o,Á COS ((a„T,. -f- 0) -1- (h) 

Rearrumando estas equações para 

sen (cu„T, + 0) -• sen 0 = 0 (i) 

COS ((W,T, + 0) - COS 0 = - 0) 

notamos que estas equações podem ser reescritas como 


sen cos 


sen sen 


+ 0 ) == 0 

(í¥+í)-s 


2(ú„mÁ 


Uma vez que sen cj„t,-/ 2 nao pode ser zero para r,- arbitrário, a Eq.^ (i’) é satisfeita 
somente se 


A Eq. (j’) torna-se então 


cos ^ + 0 j 0 
íen + 0 ) = ' 


sen ^ 

2 2(0 jn A 






de onde se tira a amplitude 

.1=_ í _ 

2a)„«isen^^ 

Pode ser de interesse o valor máximo da força da mola = kA, t nesse caso a 
Eq. (k) toma a forma 


hh 

F 


2 



Nestas condiçoes, a amplitude ou força da mola torna-se infinita quando 


co„x, 

~ 2 r 


^ = 0, n, 2n, 3n:,... 

Ji 


( 1 ) 


A Eq. (1) mostra também que a força máxima da mola é um mínimo quando 

nf„ n 3n Sn 

/; z’ 2 ’ 2 ’ ■ ■ ■ 

t\ variação de tempo do deslocamento c da velocidade pode aparecer como na 
Eig. 4.3-0. 


Pode-se aplicar um processo semelhánte, quando é -incluído 0 amortecimento, 
embora o trabalho numérico seja aumentado consideravelmente. 




Figura 4.3-6. Deslocamento e velocidade. 


4.4 FORMULAÇÃO DA TRANSFORMADA DE LAPLACE 

O método da transformada de Laplace para resolver a equação diferencial fornece 
uma solução completa, abrangendo a vibração transiente e a forçada. Apresentamos 








no Apêndice um resumo da teoria deste método destinado aos nSo familiarizados 
com 0 assunto. Nesta seção ilustraremos o seu uso com alguns exemplos simples. 

Exemplo 4.4-1 

Formular a solução da transformada de Laplace de um sistema mola-massa 
viscosamente amortecido com as condições iniciais j:( 0) e S: (0). 

Solução: A equação dc movimento para o sistema excitado por uma força arbi¬ 
trária F{t) é 

mx + cx -b kx = F(i) 

Adotando a sua transformada de Laplace, temos 

/n[í^x(.s) - x(0)í - x(0)] d- clx(s) -|- x(0)] kx{s) = F{s) 

Resolvendo em relação a x (s), obtemos a equação subsidiária 


x(s) ^ — _íl£L__ .1. ( m-r d- c)-v(O) -b mx(0) 
ms^ b CJ b á: ms^ -b CJ -1 k 


(4.4-1) 


O inverso da Eq. (4.4-1) nos dá a resposta x(r). O primeiro termo representa a vi¬ 
bração forçada, e, o segundo, a solução transiente devido às condições iniciais. 

A equação subsidiária para o caso mais geral é expressa da seguinte forma 




onde .4(í) c B{s) são polinómios, sendo que o segundo é geralmente de ordem mais 
alta que /4(s). 


Se é apenas considerada a solução forçada, podemos definir a transformada 
impedância 


Pis) 

x{s) 


z(i) 


ms^ 


cs 


k 


(4.4-3) 


A sua recíproca é a transformada admitância 

■ rij) (4.4-4) 

Usa-se freqüentemcnte um diagrama de bloco para significar a entrada e a 
saída, conforme a Fig. (4.4-1). Então a transformada admitância H{s) pode também 

Entrada F {s) - jtf (s)| -- Saída jc (s) 

Figura 4.4-1. Diagrama dc bloco. 


ser considerada como função de transferência de sistema, definida como a propor¬ 
ção no plano subsidiário da saída sobre a entrada, com todas as condições iniciais 
iguais a zero. 
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Exemplo 4.4-2 

Uma massa m é acondicionada numa caixa, conforme a Fig. 4.4-2, e cai de 
uma altura h. Deseja-se determinar a 



Figura 4.4-2. Teste de queda de uma massa acondicionada. 

a força máxima transmitida à massa meo necessário espaço de trepidação 
(rattle space*). 

Solução: Idealizamos as seguintes hipóteses; (I) Amassa m c suportada dentro da 
caixa por uma mola linear de rigidez k Ib/pol. (2) A massa da caixa é grande se 
comparada à da massa acondicionada m, de modo que a queda jivre d^ caixa não é 
influenciada pelo movimento relativo da massa m. (3) Ao tocar o piso, a caixa per¬ 
manece cm contacto com ele. 

Fazendo x o deslocamento ‘dc m cm relação à cabta, medido de cima para 
babco, a partir da posição de equilíbrio estático, e y, o deslocamento da caixa da 
posição dc partida, a equação-geral do movimento de m é 

nt(.x + .P) = —kx (4.4-5) 

Com co^ = kjm, esta equação torna-se 


Para as condições iniciais x(0), x (0), /(O), f (0). a transformada de Laplace 
para a equação acima é 

,(.)wo) I- I WO) .|. ^ <«-7) 


c a sua inversa e 


• R. D. Mindlin, “Dynamics of Packagc CushÍoning“, Bell Sy.st. Tech., Jour., 24, (julho 
1954) pa'gs. 353-461. 
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(4.4-8) 


x{t) - [x(0) -1- XO)] COS mj + ^U(0) -f- XO)] sen o)j ~ 

Podemos agora adaptar esta equação às condições do nosso problema. Dc interesse 
particular são o deslocamento -í(fo) e a velocidade x(ta') dc m no tempo to 
quando a caixa atinge o piso. 

As condições iniciais para o intervalo da queda livre são x(0) = /(O) = 
= ;i: (0) = y (0), e o movimento da caixa e a sua transformada são 

Í’(í) = ^ (4.4-9) 

Levando-so à Eq. (4.4-8) os valores da Eq. (4.4-9), encontra-se o movimento de 
m durante a queda livro, com as condições iniciais zero, o que leva a 

Sendo íq — Vs/r/? o tempo da queda da altura /i, chegamos às seguintes quanti¬ 
dades dc interesso 

•'('o) ' ■ 


aXo) -'~-scnco„r„ 

Essas quantidades tornam-se as condições iniciais para a segunda fase do problema, 
após 0 impacto da caixa sobre o piso. 


Redefinindo o tempo a partir do instante do impacto, sao as seguintes as condi¬ 
ções iniciais para a segunda fase do problema 


4(0) 

- 0, 

.v(0) - 


é(0) . 


.v(0) 

— ~ scnoj,,ro 


De acordo com a equaçao geral, Eq. (4.4-8), a equaçao para o deslocamento de 
m depois do impacto torna-se 

x(t) —^ (1 — COS CúJ) COS otj I- (--8-íq sen w | sen oj t 

ú)„ \m„ a>: " ° I 


oj; ^ 


COS co„í„y -|- {(üJa — sen scn(a)J - c^) 


(4.4-11) 


onde 



cos 0 )J„) 
-sen oj„t„) 


Assim, a amplitude máxima atingida por m é 


(4.4-12) 

(4.4-13) 


que ocorre no tcr.po (u;„ti - <t>) = 7r/2. A força máxima c simplesmente AA',. 
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Xi tem um valor máximo quando 6J„ 0, para qualquer altura do queda ■ 

h, ou tempo de quedà íq = y/2h/g. Isto se mostra recolocando sen co„fo e cos w„ío 
com suas formas em séries c admitindo que w„fo se aproxima de zero. A 
Fig. 4.4-3 mostra uma resposta de deslocamento de AT] em função da freqüência 
/„ = (l/ln) \Jkjm para /t = 10, 5, 1 c 0,15 pol.. A Eq. (4.4-13) indica, entre¬ 
tanto, que cJ^Xilg = -áfi/éjT é uma função somente de 

ct>Jo = ct>„,y2hfg ~ V2/t/<5„ (4.4-14) 

de modo que as curvas da Fig. 4.4-3 são representadas graficamente por uma tinica 
curva não-dimensional como indica a Fig. 4.4-4. 



Frcqücncia, (cps) 

Figura 4.4-3. Resposta de deslocamento para um teste de queda. 



Figura 4.4-4. Resposta de deslocamento não-dimensional 
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4.5 ESPECTRO DE RESPOSTA 


Um choque resulta dá aplicação repentina de uma força, ou outra forma de rompi¬ 
mento, que provoca uma resposta transiente de um sistema, O valor máximo da 
resposta é uma boa medida da severidade do choque e é, obviamente, dependente das 
características dinâmicas do sistema. Com o objetivo de classificar todos os tipos de 
excitações de choque, escolhe-se como um sistema-padrão um oscilador (sistema 
mola-massa) não amortecido, de um simples grau de liberdade. 

Engenhéiros acharam útil para projetos o conceito de espectro da resposta. 
Um espectro de resposta é uma representação gráfica do pico máximo da resposta 
do oscilador de um simples grau de liberdade, em função da frequência natural do 
mesmo oscilador. Diferentes tipos de excitações de choque resultarão em diferentes 
espectros de resposta. 

Considerando que o espectro de resposta é determinado a partir de um simples 
ponto na curva resposta-tempo, o qual representa um dado incompleto de informa¬ 
ção, ele sozinho não define a força do choque. De fato, é possível que espectros 
de resposta muito semelhantes corrrespondam a duas excitações de choque diferen¬ 
tes. Apesar desta limitação, o espectro de resposta é um conceito útil extensivamen¬ 
te usado. 

A resposta de um sistema à excitação arbitrária f(t) era expressa em termos do 
impulso resposta g(t) pela Eq. (4.3-1) 

= J/CÍMt - í) ifí (4.5-1) 

Para o oscilador não-amortecido de um simples grau de liberdade, temos 

de modo que a resposta de pico a ser usada no gráfico do espectro de resposta é dada 
pela equação 

■—j /(í) sen cü„(l - í)«^í (4.5-3) 

No caso em que o choque é devido a um repentino movimento do ponto de supor¬ 
te, f(t) na Eq. (4.5-3) é substituído por -ji(r), a aceleração do ponto de suporte, ou 

- — ÍXÍ)sen«„(/- (4.5-4) 

Algum tempo característico ti, tal como a duração do pulso do choque, está 
associado com a excitação do choque /(/) ou -y(t). Com r como a frequência 
natural do oscilador, o valor máximo de x(r) ou z(í) é representado graficamente 
como uma função de íj /n. 


As Figs. 4.5-1, 4.5-2 c 4.5-3 representam' os espectros de respostas para três 
diferentes excitações. A escala horizontal é 



0 0,5 1,0 Ijr 

Figura 4.5-2. 



0 1 2 3 4 5 6 

'l/r 

Figura 4,5-3. 


é igual à relação í,/r, enquanto a escala vertical é um número não-dimensional, que 
é uma medida do efeito dinâmico sobre o de uma carga estaticamente aplicada. O 
fator dinâmico de um choque é então menor que dois. 

Espectros de pseudo-resposta. Em situações de choque de solo, é muitas vezes con¬ 
veniente expressar os espectros de respostas em termos de espectros de velocidade. 
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Os espectros de deslocamento e aceleração podem então ser expressos em termos de 
espectros de velocidade, dividindo-se ou multiplicando-se por ío„. Tais resultados 
chamam-se pseudo-espectros, uma vez que eles são exatos somente quando a respos¬ 
ta máxima ocorre após haver passado o pulso do choque, em cujo caso o movimento 
é harmônico. 

Os e.spectros de velocidade são usados extensivamento cm análise de terremotos, 
e 0 amortecimento é incluído goralmentc. Com o deslocamento relativo z = x ~ y, 
é a seguinte a equação para o oscilador amortecido 

z + 2 Ííb„í -1- a>\ = -y (4.5-5) 

e a Eq. (4.5-4) é substituída por 

í j;(^)e-t“'"<'-‘:>scn «„(/ - í) (4.5-6) 

"V ^ C J 0 

Diferenciando, usando a equação 

Jj-1' /(Í.í) r/í - f r/í 1- /(í, í) í í (4.5-7) 

obtemos para a velocidade 

[—ÍCü„ sen v'! —X^ ®„(/ — í) 

+ <»„V1 — COS VI - í‘ w„(r — í)] r/í (4.5-8) 

Fazendo 

-4 = í COS VI — í^ «„í r/í (4.5-9) 

J 0 

(4.5-10) 


(4.5-11) 
ou 


i? = í' Kí/e''""' sen V w„í r/í 

J 0 

A Eq. (4.5-8) pode ser expressa como 

i(0 = ; y i[[-4; - 5yn=y^]scnVr^'^ra„r 
+ [/ivr=y -I- BQ COS yr^ ca„í} 



Se a Eq. (4.5-12) é traçada cm função do tempo, ela aparece como uma onda 
modulada de amplitude, como se vê na Fig. 4.5-4. Assim, a velocidade máxima de 
resposta S^,, ou o espectro de velocidade, c dada com exatidão suficiente pelo valor 
pico da envolvente 


Figura 4,5-4. 





Para a segunda função-rampa começando em pelo exame da equação aci¬ 
ma, a solução pode ser expressa como 

x(,\- senm,(/ -r,) T 

“ kl t, “ J 

Superpondo estas duas equações, a resposta para t > torna-sc 


x{l) = ^[l 


sen cü / ,• I , „ 

.-2_ J-senct; 

CúJi CúJi 


-1,) 


Diferenciando e igualando a zero, obtém-se o tempo pico como 

tg 

^ " " sen co^í, 

Uma vez que ío„tp deve ser maior que rr, obtemos também 
sen coj^ ~ —1 — COS coj^) 

—sen CO t. 

COS = - ■ ” ‘ 

V2(l — COS cü„í,) 

Substituindo essas quantidades em x(í), acha-se a amplitude pico 


if) - 


/2(1 — COS 


Fazendo-se o período do oscilador r = 27r/cj„, a equação acima é traçada como 
função de ti /r, como na Fig. 4.5-1. 


Exemplo 4.5-2 

Determinar o espectro de resposta para a entrada da velocidade base, y (t) = 
= Vç,e~'^‘i do Exemplo 4.3-2. 


Solução: O deslocamento relativo z(/) foi achado no Exemplo 4.3-2, .sendo 
^(0 = 1 ^ '* ~ "-'o^en coj - COS coj) 

0 processo usual para determinar o valor pico de é diferenciar a equação em re¬ 
lação ao tempo f, igualá-la a zero, e substituir este tempo dc volta na equação 
para z(í). É evidente que a conseqücncra neste problema é uma equação transcen¬ 
dental que tem de ser resolvida graficamente. Para evitar esta tarefa numérica, vamos 
considerar, a seguir, um processo diferente. 

Para um sistema muito rígido, que corresponde a grande, a resposta pico 
ocorrerá certamente para t pequeno, e obteríamos para a parte da equação que varia 
com 0 tempo o valor pico de 

(1 - mja - 1 ) = -cojo 

Deste modo, para grande, o valor pico será aproximadamente igual a 


EàLs- 


de modo que o traçado de Zplvoíg versus é um hipérbole retangular. 


Figura 4.5-6. Doblele impulsivo. 


Para co„ pequeno, ou uma mola muito flexível, a duração da entrada seria pequena 
comparada ao período do sistema. Daí, a entrada apareceria como um doblete im¬ 
pulsivo, conforme a Fig. 4.5-6, com a equação Voto^XO- Então, a solução para 
z(f)é 

z(l) — rJoto COS Cüj 

e o seu valor pico é 

UpISVo 

Avaliadas estas condições extremas, podemos agora preencher o espectro de 
resposta, que é indicado na Fig. 4.5-7. 



4.6 O COMPUTADOR ANALÓGICO 

O nosso breve encontro com tránsientes e cspectrqs de resposta mostra suficiente¬ 
mente as dificuldades algébricas com que deparamos, .mesmo nos problemas muito 
simples. Tais problemas são resolvidos de modo adequado pelo computador analó¬ 
gico, que serve idealmente para a solução de equações diferenciais comuns. Esses 
computadores são capazes também de resolver problemas não-lineares, para os quais 
as técnicas analíticas são inexistentes ou complexas demais para uso prático. 
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0 elemento básico do computador analógico é o amplificador operacional de 
■alto ganho, indicado simbolicamente na Fig. 4.6-1. Não c, entretanto, necessário 
conhecer os detalhes da sua parte eletrônica. Tais amplificadores são caracterizados 
pela equação 



Figura 4.6-1. Amplificador operacional. 

onde fi é Q fator de amplificação e e' Cq são as voltagens da rede e da saída. O 
fator de amplificação /r para os amplificadores operacionais modernos é apro¬ 
ximadamente 10®. Considerando que a voltagem de saída c limitada geralmentc a 
± 100 volts, a ordem dc magnitude da voltagem da rede c ± 10'® volts. A cor¬ 
rente consumida pela rede é também muito pequena, com um valor representativo 
de 10'^ amps. 

Várias operações de diferenciação, integração, soma, etc., são possíveis dc 
realizar, pela conexão do amplificador operacional a diferentes tipos dc impedância. 
A Fig. 4.6-2 indica um entrelaçamento geral do amplificador com uma impedância 



^ Figura 4.6-2. Circuito do amplificador operacional. 

) de entrada Zi e uma impedância de rcalimentação Zj-. As seguintes equações são 
^ relativas ao circuito acima 

) e, — e, = í,Z, 


="= ‘e -I- í. 


‘'o 


Nestas equações, é insignificante'em comparação com Co e e- c í é uma 
quantidade desprezível em comparação com q e U Assim, com e”= i d as 
equações acima sc transformam em s g ’ 


-í’o — i/Z^ 


Obtemos, a partir delas, a relação entrada-saída 

<^0 Zf 

(4.6-2) 

que c fundamental para o computador analógico. 

(n) Mudança dc sinal. A .nudança de sinal é a operação mais simples. Fazen¬ 
do z,. = R. c Zf = Rj., a Eq. (4.6-2) torna-se 

A, 

(4.6-3) 

de modo que se Ay = R. obtemos simplesmente 

Co — —Cl 

A Fig. 4.6-3 mostra o circuito para a mudança de sinal. 



Figura 4.6-3. Circuito para o fator dc escala R^/R-, 

..do ““““ 

cado na Fig. 4.6-4, então /y e a soma das correntes de entrada 

h d" '2 “)■ I3 -- d 

^ -1- . e„ 

■R. /<2, ' Rt 


e a voltagem dc saída é dada pelo somatório 

p . _V 

” ' . (4.6.4) 

Se mdas as resistências são iguais, a consequência é o somatório das voltagens de en- 
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Figura 4.6‘4. Circuito para soma. 


(c) Integração. Se a impedância de realimentação c um capacitor C, como 
indica a Fig. 4.6-5, o circuito efetuará a função de integração. Com a voltagem 
inicial e(0) através do capacitor, a sua voltagem a qualquer tempo t é (lembran¬ 
do que s 0) 



Figura 4.6-5. Circuito de integração com condições iniciais. 

Porém, i = e,/K, de modo que a equação acima torna-se 

<•„ =.- '7^ [ ^ <-(0) , (4.6-5) 

Com R = 1 megaohm c C = 1 microfarad, RC = 1 segundo, e o tempo do 
computador é diretamente em segundos. 

A voltagem inicial e(0) c obtida do circuito indicado em linha pontilhada, 
fechando-se a chave S antes do iniciar a computação. Iniciada esta, a chave S 6 
simultaneamente aberta por um relé. 
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(d) Diferenciação. Não se faz diferenciação cm computadores analógicos pelo 
fato do inevitável sinal de ruído na entrada ser ampliado por p, provocando a satu¬ 
ração do amplificador. Em vez, é usualmente possível rearrumar as equações para 
integração. 

(e) Divisão de voltagem. O potenciômetro muitas vezes usado para se obter 
uma fração k vezes a voltagem de entrada 


é designado simbolicamente na Fig. 4.6-6 (a). O ajustamento fracionário k no po- 
tenciômetro só prevalece quando a saída está em circuito aberto. 




Figura 4.6-6. Divisão de voltagem. 

Quando uma resistência de carga R^ c colocada através a saída, como na Fig. 4.6-6(b), 
pode-se mostrar que a voltagem do saída é igual a 


ke, 


1 -I- ^ki\ - k) 


onde A é a resistência do potenciômetro. É evidente que esta equação se aproxima 
de kci quando RjR^ 0. 

(/) Multiplicação. Multiplicação é uma das mais difíceis operações para o 
computador analógico. Em um método, o princípio do potenciômetro descarregado 
é aproveitado usando-sc um scrvomccanismo com um potenciômetro conjugado, 
conforme a Fig. 4.6-7. O primeiro potenciômetro c ligado a ± 100 volts, ao passo 
que o segundo c ligado a te^, a voltagem a ser multiplicada por Cj. 

+ 100 +«! 


e, “ lÜOA: Servo- 

-^—>- ir)cca* 

nismo 


Figura 4.6-7. Servomultiplicador. 










1 


A função do servomccanismo é posicionar - o cursor do potcnciômctro cali¬ 
brado - em zero o erro entre a saída do primeiro potcnciòmetro c a voltagem de 
entrada «i. Visto que a saída de cada potenciômetro é proporcional a k, temos 

(7, = lOOA- 0 Co — Avj 

Com a eliminação de k, obtemos a equação 


ig) Circuito do Computador ^ara o Sistema de Um Grau de Liberdade. O cir¬ 
cuito da Fig. 4.6-8 demonstra o uso do computador analógico na solução do sistema 
linear de um grau de liberdade. A ecíuação representada c 

m.K - c.x -[- k.x ■— Lil) 


1 


'■5, 

x(0) ç 

-X R 

^'c~' 

LKi 


(ici) = 2 f u 



Figura,4.6-8. Circuito anaiógico para p sistema de um grau de liberdade. 
que é rearrumada para 

.V == —2fcü„Á' — Cúl-x -f -^^(0 

Admitindà que a entrada do primeiro amplificador seja x, sua saída i —x, 
etc. As voltagens nos três terminais © , (5) , c @ , que são iguais a x são aque¬ 
las do lado direito da equação acima. Observar que o sinal muda através dc cada ampli¬ 
ficador e que os ajustamentos do; potenciômetro são para 2fco„ e oj^. As condi¬ 
ções iniciais x(0) e ,x:(0) são as voltagens nos dois capacitores no tempo r = 0. 
Quando o computador é posto em funcionamento pelo fechamento da chave S^, as 
chaves t S-í abrem-se simultaneamente. 

Qi) Mudança de Escala. Sc as variáveis do problema passam por urna trans¬ 
formação linear, as características do sistema permanecem as mesmas, e apenas ex¬ 
pandimos ou contraímos a escala das variáveis. Na solução de um problema no 
computador, é muitas vezes necessário fazer tal mudança de escala, em razão das li¬ 
mitações do computador. Uma mudança na escala do tempo pode ser necessária pelo 


fato da vibração efetiva poder ser de freqüência alta demais para o computador e o 
registrador acompanharem. Por outro lado, se a vibração efetiva c baixa demafs cm 
freqüência, o tempo do computador torna-se longo cm excesso, e por isso introduzin¬ 
do erros por desvios. O escalamcnto da amplitude é também necessário, de modo que 
os amplificadores operem dentro dos limites de ±100 volts. A resposta pico deve 
ser aproximada dos limites dc ± lOO-volt, para o máximo de exatidão. 

Suponhamos que a equação para o sistema vigente seja 

-v(r) I- 2ííu„Ã-(r) -I- (o^xii) (4.6-8) 

com as condições iniciais 

x(0) e A-(0) 

Sc desejamos niudar dc a a escala de tempo do problema, fazemos 


Então as equações seguintes dão as derivadas „ 


dr^ dE 


(Ít iít dx^ dt^ 

e a equação diferencial Original com.suas condições iniciais torna-se 

-|. colxix) : = If© " ' 


1 r/.v(0) , 
a dt'' 


.v,(0) .v,(0)' 


(4.6-9) 


(4.6-10) 


(4.6-11,) 


Dividindo-se por , esta equação toma a forma dc ■. ‘ 

que mostra que a freqüência natural do sistema mudou de para Cl = co„/a. Ó 


fator de amortecimento f, entretanto, não mudou, uma vez que o amortecimento 
crítico = 2m_a)„ para a equação original mudou para = ImCl para a 
equação nova. 

Teoricamente, c possível planejar um circuito de computador para resolver 
esta nova equação c interpretar seus resultados em termos das variáveis priginais. 
Entretanto, há problemas relativos às ordens de magnitudes das voltagens do circuito 
do computador que precisam dc mais atenção, c cuja discussão ê melhor em termos 
do exemplo seguinte. ' \ \ 

Exemplo 4.6-1 > 

A equação para certo sistema mecânico, excitado por urna carga degrau /(í) = 
= 2000 Ib c apresentada como ■ 

0,10.v. 5.V -f 4000x = 2000 Ib, . 

com as condições iniciais 





^■(0) = -20pol/s e .x(0) = 0.25 pol. 

Hscrever novanientc a equação para o computador e estabelecer um circuito 
praticável para a sua computação. 

Solução: Escrevendo a equação na forma 

.v + 50i + 40,000a: = 20,000 
encontramos para a freqüéncia natural do sistema o valor 

(ú„ = ^kjm = ^"40,000 = 200 rad/s 

Esta freqüéncia é muito álta para o computador, de modo que escolhemos arbitra¬ 
riamente a = 100 para baixar por um fator de 100 a escala do tempo. A nova 
equação com t como variável independente é então [Vide Eq. (4.6-12)] 

-1- 0,50?^ -b 4,0 .y(t) - 2,0 
dx 

com as condições iniciais 

^ - -0.20. .v.(0) 0,25 

e observamos que a freqüéncia natural ficou reduzida a íí = 2 rad/s. 

Se ignoramos as ordens de magnitudes, o seguinte circuito da Fig. 4.6-9 satisfaz 
essas equações. Todavia, para que o computador dê resultados seguros, as voltagens 
de saída dos amplificadores não devem exceder ±100 volts, nem serem pequenas 
demais. Em face dessas considerações, é necessário fazer uma estimativa do deslo¬ 
camento, velocidade e aceleração máximos a serem encontrados na equação dife¬ 
rencial modificada e estabelecer fatores apropriados de escala, que darão voltagens 
pico de saída de amplificadores aproximadas do máximo admissível de ± 100 volts. 



Figura 4.6-9. Circuito analógico para um sistema de um grau de liberdade. 
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Através da equação modificada, reconhecemos que o amortecimento elimina¬ 
rá eventualmenle todas as oscilações e que o deslocamento final atingido será igual a 

Notamos também que, sem o amortecimento, a amplitude pico sob a excitação de 
uma função degrau é duas vezes o valor acima, ou 

x(t) = 1,0 pol 

mtx 

A velocidade e aceleração máximas a serem encontradas são as seguintes, estimadas 
na base de movimento harmônico 

=Qx(t) = 2 pol/s 

}_ dx Jm»x mix 

rr/bv(T)~j ^ Cl^xir) = 4- po!/s^ 

[_ dx Jtnix mix 

Se examinamos agora o circuito da Fig. 4.6-9 com os valores máximos acima, consta¬ 
tamos que os amplificadores 1, 2 e 3 têm voltagens pico de saída de 4, 2 e 4 volts 
somente, e que não podemos obter exatidão com tal circuito. 

Para superar esta dificuldade, efetuamos uma mudança de escala de modo que, 
com x(t) = 1 pol, a salda estará mais aproximada do pico admissível de 100 volts. 

A fim d^evitar que ultrapasse o limite 100-volt, fixaremos esta voltagem em 80 volts. 
Multiplicando a equação diferencial modificada por 20, obterhos 

-10'-^:^ - 80a-(í) 40 

í/t dx 

O circuito para esta equação de escala acrescida toma agora a forma indicada na 
Fig. 4.6-10, onde as voltagens iniciais são escaladas de acordo com a saída dos ampli¬ 
ficadores, isto é, com c/x{0)/dT = - 0,20 pol/s, e a saída do amplificador 2 ig^ial a 



Figura 4.^-10. 











) 

) 


— 10(íiç(T)/dT), seu valor seria — 10(—0,20) = + 2 volts. De .forma semelhante, a ; 

^ voltagem inicial do amplificador 3 seria 80 x 0,25 = 20 volts. ■ 

) Temos agora uma saída de 80 volts para o amplificador 3, correspondendo à ; 

amplitude máxima esperada de 1 pol. Entretanto, uma vez que a velocidade máxima I 

^ esperada é de 2 pol/s, a saída pico do amplificador 2 é somente 10(cix/dT) = 20 volts. j 

) É preferível, então, mudar o ganho do amplificador 2 por um fator de 4 e reduzir j 

consequentemente o-ganho do amplificador 3, conforme mostra a Fig. 4.6-11. Sendo j 

) de 4 pol/s o valor máximo esperado da aceleração, não há necessidade de mudar o í 

•j ganho do amplificador 1. O ganho do amplificador c determinado pela quantidade 

RC, que é unidade para- R = 1 megaolun o C = 1 microfarad. Observe-se também í 

) que essas mudanças requerem mais uma mudança na voltagem inicial do amplificador > 

X 2, dc 2 para 8 volts. ? 

Finalmonte, esses circuitos do computador aparentam não serem únicos, c a í 

^ mesma equação pode scr resolvida por diferentes circuitos. Assim, a equação para r 

i este problema pode também sor resolvida pelo circuito da Fig. 4.6-1 2. ) 




4.7 DIFERENÇAS FINITAS EM COMPUTAÇÃO DIGITAL 

Quando a equação diferencial não pode ser integrada na forma fechada, métodos nu¬ 
méricos devem ser empregados. Este pode bem ser o caso quando um sistema mo- 
la-massa é excitado por uma força que não pode ser expressa por funções analíticas 
simples, ou quando uma tal força é dada gráfica ou numericamente apenas. 

A integração numérica é um processo pelo qual a equação diferencial dc mo¬ 
vimento é resolvida progressivamente cm incrementos dc tempo, partindo dc algum 
tempo, quando o deslocamento c a velocidade eram conhecidos. A solução ó aproxi¬ 
mada mas, a medida que o intervalo de tempo é reduzido, o resultado se aproxima 
da solução exata. Embora existam diferentes métodos numéricos disponíveis, vamos 
considerar neste capítulo apenas dois, escolhidos pela sua simplicidade. Os méritos 
dos vários métodos são associados aos erros e convergência do processo. Estes são 
discutidos cm muitos trabalhos sobre análise numérica.* 

A base da solução numérica c csscncialmente.a de se obter os valores numéricos 
da integral, a scr determinada, nos pontos pivotais ao longo do eixo do tempo. Neste 
sentido, as derivadas na equação diferencial são aproximadas dc um certo número de 
termos na expansão dc Taylor. Na expansíp Taylor, ãr,-.,., e xr,- ,, cm termos 
dos pontos pivotais são 


-V,-M 

.V,. 

i -V, At 1 





- .V, A< 1 




Subtraindo e não considerando os termos dc ordens mais altas, obtemos 
. .f.Ar .v,.,) />;2 

Somando, achamos 


Estas equações recorrentes, juntamente com as equações diferenciais de movi- 
mento,, são suficientes para a solução numérica: Entretanto, algumas considerações 
sao necessárias para começarmos o processo de computação. 


A. Kal.ston c II.‘S. Wiif, Malhcnmtical Mcthods for Digital Computer,'; Vols I & II (Nova 
Iorque;iohn Wiley & Sons, 1968), , , ■ : 

1952) Sal.vadori c M. L, Baron, Numerical Mcthods, in. Fngincçring (Prcnticc-Hall, -Inc., 

Lll. 








Quando a aceleração inicial (ou força) não é zero, o processo mais simples é 
admitir que ela permaneça constante durante o primeiro intervalo. Como subzcro 
não é disponível no computador digital, usamos sub 1 para os valores iniciais. Temos 
então 


.Vj-=y.V,Aí^ 

A aceleração x, a ser usada na equação acima é determinada a partir da equação 
diferencial de movimento e suas condições iniciais Xi e x i 



O deslocamento X 2 e velocidade x^ são então detenninados pela Eq. (4.7-4) 
e substituídos na equação diferenciai para üj. Tanto x^ c x^ são então substi¬ 
tuídos na Eq. (4.7-3) para determinar .Xj e o processo se repete. Usualmente, 
um intervalo de tempo At < l/IOr, onde t é o período natural, é suficieniemente 
pequeno para resultar numa solução satisfatória. 


Um diagrama de fluxo para o cálculo digital é indicado na Fig. 4.7-1. 



figura 4.7'I. 


Cora os dados fornecidos no bloco ® seguimos para o bloco ® que é a equação 
diferencial. Indo a © pela primeira vez, / não é maior que I, e daí seguimos para 
a- esquerda onde Xj é calculado.. Aumentando /'de I completamos a iteração 
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esquerda para (g) e © onde 1 é agora igual a 2, e assim seguimos para a direita 
para calcular Xj. Admitindo N intervalos de , Ar, o caminho é para a direção 
NO e a iteração direita é repetida N vezes até / = A' + 1, cm cujo tempo os 
resultados são impressos. 

O programa Fortran para a computação é dado conforme 0 diagrama de fluxo. 


FORTRAN PROGRAM 


/= 1 


20 

X, = f{x„ 1 

+ (/ - 1) Ar) 


IF (I.CT.I) 

GO TO 30 


.Xj |.v, At 

^ + x, At -1- X, 


/ = / -M 

GO TO 20 

30 

X, == 2.x,,, 

1 

1 

> 


IF (/ =- N 

I - I -1- 1 

GO TO 20 

• 1) GO TO 40 

40 

PRINT 



STOP 



END 



Se a aceleração inicial xi c zero, a Eq. (4.74) resulta em Xj = 0, e 0 pro¬ 
cesso de cálculo não pode corneçar. Esta condição pode ser retificada pelo desenvol¬ 
vimento de uma equação nova baseada na hipótese de que a aceleração varia linerar- 
mente de Xj para Xj durante o primeiro intervalo, como se segue 

.V, - = .Vr -1- a Aí 

.X,- ■ .V, Ar i ya Ar- 

.V, ‘^-^.v,Aí= -l--^aAf’ 

Eliminando a na última equação, tirando o seu valor da primeira, obtemos 


(4.7-5) ■ 


-^-(2.V| i .V,) A/= 


Com Xi = 0, a Eq. (4.7-5) e a equação diferencial têm de set resolvidas por tenta- 


tivní í» • 
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(4.7-6) 


1 

> 

. ' 

, > 

^ ■ X2 <!>(■'■'1’ ^2) 

• ^ 0 método Runge-Kutta, que é auto-iniciado, é discutido na Seção 4.8. 

I ) 

Exemplo 4.7-1 

I Resolver numericamente a equação diferencial 4x + 2000a' — F{í), com 

j ) as condições iniciais Al = Á, = 0 c l‘{t) como estabelecido na lúg. 4,/-2. 

I j 

- 100 

I ) 

' ) 

' ) 0 ■■ 0,10 0,20 t 

I ) Figura 4.7-2. 

^ Solução: Encontra-sc primeiro o período natural do sistema, que é 

i) ^ ' cu==^- 7^.= 22,4 rad/s 

' ) 0.281 s , 

I Neste problema utilizaremos a régua de cálculò para efetuar a computação, e esco- 

' lhemos t = 0,030, que é aproximaflamentc 1/10 t, a fim de reduzir a um mínimo 

* ^ 0 número de computações. .. 

!• . A equação diferencial nos dá a aceleração qüe é 

' . ' ' ■ ■ -é - 500.V ' , . , 

) 

' e obtemos, com Xj - 0 

I ^ -í, = i.x 100 - 0 ^ 25 

I ^ A Eq. (4.7-4) dá então 

, I Al -1(25 X0,030)^ =0,0113 

i ) Com A, = 0,0113 voltamos à equação diferencial para determinar. A, . 

) , .Xi =-1(100) r 500(0,0113) 

■ = 25 - '5,65 = 19,35, 

. As quantidades x, e são agora substituídas na Eq. (4.7-3) com / = 2 
' ^ ' 19,35(0,030)^ --= 0 - 2(0,0113)-r Aj 

' ■ =i0;0400' ' ' ■ ' ' . , \ 

) 1,14. 

) 



) 


Temos agora x, e x, ; x, pode ser computado a partir da equação diferencial. 
A Eq. (4.7-3) dá então x^, c o processo é repetido entre a equação diferencial e a 
Eq. (4.7-3). A tabela seguinte ilustra a sequência do cálculo, cujos resultados são 
representados graficamente na Fig. 4.7-3. 


Figura 4.7‘J, 


/ 

t 

VV) 

5Ü0.V 

X 

Xàt~ 

.V 

1 

0 

25,0 

0 

25 

0,0225 

0 

2 

0,03 

25,0 

5,65 

19,35 

0,0174 

0,0133* 

3 

0,06 

25,0 

20,00 

5,00 

0,0045 

0,0400” 

4 

0,09 

25,0 

36,60 

-11,60 

-0,0104 

0,0732 

5 

0,120 

20,0 

48,0 

-28,0 

-0,0252 

0,0960 

6 

0,150 

12,5 * 

46,80 

-34,3 

-0,0309 

0,0936 

.7 

0,180 

5,0 

30,15 

-25,15 

-0,0226 

0,0603 

8 

0,210 

0 

2,20 

- 2,20 

-0,0020 

0,0044 

9 

0,240 

0 

-26,75 

26,75 

0,0241 

-0,0535 

10 

0,270 

0 

-43,65 

43,65 

0,0393 

-0,0873 

11 

0,300 

0 

-40,90 

40,90 

0,0368 

-0,0818 

12 

0,330 

0 

-19,75 

19,75 

0,0178 

-0,0395 

13 

0,360 

0 

10,30 

-10,30 

-0,0093 

0,0206 

14 

0,390 

0 

35,70 

-'35,70 

-0,0321 

0,0714 

15 

0,420 

0 

45,05 

-45,05 

-0,0405 

0,0901 

16 

0,450 

0 




0,0683 

DaEq. (4.7-4) A, 

= ii.Ar‘ = 

-j(0,0225) = 

= 0,0113 



Da Eq. (4.7-3) Xj 

= -A', + 2a:, 

+ Aj Aí’ = - 

-0 + 2(0,0113) + 

0,0174 = 0,040 



Exemplo 4.7-2 

Resolver pelo computador digital o problema de um sistema mola-massa ex¬ 
citado por um pulso triangular. São as seguintes a equação diferencial e as con¬ 
dições iniciais 
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0,5.í + 8;r^ji: F{i) 
.V, = .Y, =: 0 

A força triangular é definida na Fig, 4.7-4 



Solução: O período natural do sistema é 

T ^ = 2?r 0,50 

cu 4n 

■ O incremento de tempo será escolhido como t - 0,05, e a equação diferencial é 
reorganizada como 

.V ■— 2F(!] lÓTt^.V 

Esta equação é para ser resolvida juntamente com a equação recorrente (Eq. 4.7-3) 
.v,,| = XiAF ■ - .Y,., -i- 2.xI 

A força e a aceleração sendo zero quando t = 0, é necessário começar o 
processo computacional com as Eqs. (4.7-6), que são 

Xj - i;f2(0,05)^ = 0,000417x2 

Xj = 2F(0,05) - lóa-^Xj = 50 - ISSxj 

A sua solução simultânea conduz a 

(0,05)^ F(0,05) 

^ ( = 0,0195 

3 -t 87r^(0,05)^ 

.Vj -,46,91 

O diagrama do fluxo para a computação está representado na Fig. 4.7-5. Cora 
At = 0,05, a duração do tempo para a força deve ser dividida em regiões / = 1 
até 5, / = 6 até 9 e / > 9. O índice / controla o caminho da computação no 
diagrama. 

O programa Fortran pode ser escrito de muitos modos, um dos quais está 
indicado na Fig. 4.7-6, e os resultados, na Fig. 4.7-7, podem ser represeptados grafi¬ 
camente pelo computador, conforme a Fig. 4.7-8; Um At menor resultaria em 
gráfico menos acentuado. 

16 




At = 0.05 



■MÉIi 

















c • ■ 

C VIRRATION PROBLEM 
C 

I5N 0002 DIMENSIOfl X(25},0X2(25).P(25),T{25),J(25),VAR(25) 

ISfl 0003 PI20.1415**2 

ISfl 0004 0T=0.O5 

ISfJ 0005 0T2=0T**2 

ISN 0006 X(1)=0.0 

ISÍI 0007 0X2n)=0.0 

is;i 0000 F(i)=o.o . ■ 

ISN 0009 T(1)=0.0 

ISN 0010 J(1 ) = 1 

ISN oon 00 1 1 = 2,25 

!SN 0012 J(I)=I 

rsri 0013 T(i)=OT*(i-i) 

ISN 0014 IF (I .GT. 2) GO TO 2 

ISN 0016 F(I)=500*0T*(I-1) 

ISIi 0017 ■X(I) = (0T2*F{1) )/{3+8*PI2*0T2) 

ISN 0018- DX2(I)=2*F(I)-16*PI2*X(I} 

ISN 0019 GO TO 1 

ISN 0020 2 IF(I .LE. 5) F(I)=500*DT*(I-1) 

ISN 0022- IF (I .GT. 5 .ANO. I .LT. 9) F(I)=200-500*DT*(1-1) 

ISN 0024 IF (I .GF. 9) F{I)=0.0 

ISN 0026 X{I)=ÜX2(.-l)*DT2-X(I-2)+2*X{I-l) 

ISN 0027 DX2(I)*2*F(1)-16*PI2*X(Í) 

ISN 0028 1 CONTINUE 

ISN 0029 '.;RITE(6,3) 

ISN 0030 3 F0RMAT(4TH1 J TIME OISPL ACCLRTN FORCE) 

ISN 0031 WRITE(6,4} (J(IJ,T(I).X(I),DX2(I),F(I).I=1,25) 

ISN 0032 4 F0RHAT{3X.I2,2X,F6.4,3X.F6.3.3X,F7.2.4X,F7.2) 

C 

C PLOTTINC. 

c 

ISN 0033 DO 5 1=1,25 

ISN 0034 5 VAR(I)=X(r)*10 

ISN 0035 CALl. PLOTl ( VAR ,25) 

ISN 0036 STOP 

ISN Q037 CND 

Figura 4.7-6. 


J 

TIME 

OISPL 

ACCLRTN 

FORCE 

1 

0.0 

0.0 

0.0 

0.0 

2 

0.0500 

0,020 

46.91 

25.00 

3 

0.1000 

0,156 

75.31 

50.00 

4 

0.1500 

0.481 

73.97 

75.00 

5 

0.2000 

0.992 

43.44 

100.00 

6 

0.2500 

1.610 

-104.25 

,75.00 

7 

0.3000 

1.968 

-210.78 

50.00 

8 

0.3500 

1.799 

-234.10 

25.00 

9 

0.4000 

1.045 

-165.01 

0.00 

10 

0.4500 

-0.122 

19.22 

0.0 

11 

0.5000 

-1.240 

195.86 

0.0 

12 

0.5500 

-1.869 

295.19 

0.0 

13 

0.6000 

-1.760 

277.98 

0.0 

14 

0.6500 

-0.957 

151.04 

0.0 

15 

0.7000 

0.225 

-35.52 

0.0 

16 

0.7500 

1.318 

-208.06 

0.0 

17 

0.8000 

1.890 

-298.47 

0.0 

18 

0.8500 

1.717 

-271.05 

0.0 

19 

0.9000 

0,865 

-136.64 

0.0 

20 

0.9500 

-0.328 

51.72 

0.0 

21 

1.0000 

-1.391 

219,66 

0.0 

22 

1.0500 

-1.906 

300.89 

0.0 

23 

1.1000 

-1.668 

263.33 

0.0 

24 

1.1500 

-0.772 

121.83 

O.Ü 

25 

1.2000 

0.429 

-67.77 

0.0 


Figura 4.7'7. 


F X 
100 - 2 


50h Ih 


qI- 0 x-x- 


0,2 0,4 J, -0,6 / 0,8 \l,0 1,2/ 



Figura 4.7-8. 


4.8 A COMPUTAÇÃO RUNGE-KUTTA 

0 processo de computação Runge-Kutta é muito popular, por ser auto-iniciado e 
apresentar resultados com precisão satisfatória. A seguir, uma discussão resumida 
da sua base. 

Consideremos a equação diferenciai para o sistema de um grau de liberdade 

iPx , d.X , , ’ rr.\ 

Fazendo / = dx/dt, a equação acima pode scr escrita como 

(4.8-1) 

X e y, ambos na vizinhança de Xj e podem ser expressos em termos da série 
dcTaylor. Fazendo o incremento de tempo h = At 



m 




Em lugar de usar estas expressões, é possível substituir a primeira derivada por 
uma inclinação média e ignorar as derivadas de ordem mais elevada 














.V ... ,x, -1 (^) h (4.8-2) 

V -.v/ l (4.8-3) 

Se tivermos usado a regra de Simpson, a inclinação média no intervalo h torna-se, 
isto é 



O método Runge-Kutta é muito semelhante às computações acima, com a 
exceção de que o termo central da Eq. (4.8-4) é dividido em dois termos, e quatro 
valores de !, x, y, e. f são computados para cada ponto como se segue 

, .V V .t f ’ y - x 


r, /, 


X, - .V,. 


V-, ,iv 


F, - 

’fO'uX\. y,) 

Ti - í( 

* 2 

-V; - .Vi 

, y h 
* 1 ^ 

y, .... . 

. rr 
i Tiy 

Fi - 

■ f{Ti, Xi. Yi) 

Ti 1, 

1-- 

,Vj .V, ■ 

, y 

t- Yi-y 

>’j - .'V 

. /.' 

r- / 2 — 

Fi - 

fiTi,Xi, Yi) 

Ti - 1, 

1 /l 

.V4 ■ .V, 

1 Y,li 

)'4 ^ - .Vi 

i Tih 

F, 

J(Ti, Xi. Yi) 


Essas quantidades são então usadas nas seguintes fórmulas do recorrência 
. .v,.M ■ 1 4 ^^' ■' 

V,., ...v, .| 4[F, -1 2F, -1 2F, (4.8-6) 

onde SC reconhece que os quatro valores de Y divididos por 6 representam uma in¬ 
clinação média dxjdt e os quatro valores de F divididos por 6 resultam numa 
média de dxjdt, como definidas pelas Eqs. (4.8-2) e (4,8-3). 


PROBLEMAS 


4-1 Mostrar que o tempo í^, correspondente à resposta pico para o sistema 
mola-massa excitado impulsivamentc, c dado pela equação 


4-2 


tg \/1 - cu„t, - - V1 — cvc 

Determinar o deslocamento pico para o sistema mola-massa impulsivamentc 
excitado, e mostrar que ele pode ser expresso na forma 

Representar graficamente este resultado como uma função de f. 


1 4-3 Mostrar que o tempo íp correspondente à resposta pico do sistema mola-n^- 


sa amortecido, excitado por uma força degrau I'a, é co„tp 
4-4 Mostrar que, para o sistema do l’robl. 3, a resposta pico é igual a 

írt 


rr/Vr:?. 


(â. 


4-5 Um pulso retangular de altura /'o e duração íq é aplicado a um sistema 
mola-massa não-amortccido, Considerando ser o pulso a soma de dois pulsos 
degraus, como indicado na Fig. P.4-5, deternunar sua resposta para t > to 
pela superposição das soíuçocs não-amortccidas. 


Figura 

4-6 Sc uma força arbitrária f{t) c aplicada a um oscilador não-amortccido, que 
tem condições iniciais diferentes de zero, mostrar que a solução deve ser na 
forma 


.V(,) :..v„cosw„t I ^senov I 


/(í)senoa„(í — í) t/í 


4-7 Mostrar que a resposta a uma função degrau unitária, designada por /i(r), 

c relacionada à resposta impulsiva g(l) pela equação g{t) = /i(í). 

4-8 Mostrar que a integral de convolução pode também ser escrita em termos de 
li(t) como 

.V(í) f(OViU) i I f(OliU-OdÍ 
0 

onde h{t) c a resposta a uma função degrau unitária. 

4-9 Na Seção 4.4, a Eq. (4.4-1) dá a equação subsidiária para o sistema mola-massa 
viscosamente amortecido. Avaliar o segundo termo devido a condições imciais 
pelas transformadas inversas. 

4-10. Um sistema mola-massa não-ámortecido recebe i^ma excitação base de j>(/) = 
= 90(1 — 5/). Determinar o deslocamento relativo máximo, sabendo-sc que 
a freqüência natural do sistema é = lOs . 

4-1! Um pulso scnoidal é considerado como a superposição de duas ondas senoidais, 
conforme a Fig. P.4-11. Mostrar que a solução é . ^ ^ 


20 



) 

) 



Figura P.4'U. 


4-12 Mostrar que a resposta para o pulso triangular indicado na Fig. P.4-Í2 é 



—2-(f-l/2í,) 

í, 

Figura P-4-t2. 


^ - T'') - ,,) -scn2;r-f]}, / > ,, 

4-13 Um sistema mola-massa desliza para baixo num plano inclinado suave de 30°, 
conforme a Fig. P.4-13. Determinar o tempo decorrido entre o primeiro con¬ 
tato da mola até que ela desfaça novamente o contato. 



Figura P.4-13. 


4-14 Com referência ao Exemplo 4.4-2, determinar o espaço de trepidação (“rattle 
spacc”) necessário, no ca.so do .sistema suspenso ter uma frequência natural 
dc 10 cps c da caixa cair de uma altura de 1,0 pól. 

4-15 Um peso dc 38,6 Ib c suportado por várias molas, cuja rigidez combinada é 
dc 6,40 Ib/pol. & o sistema c suspenso de modo que a base dar. molas é ao mes¬ 
mo tempo livre c solta, detenninar o deslocamento máximo dc m, o o tempo 
para a compressão máxima. 

4-lc Um instrumento delicado é suportado por molas dentro de uma caixa, confor¬ 
me indicado na Fig. í’.4-16, onde sua frequência natural é de 10 cps. Na posi¬ 
ção dc suportado, o espaço livre entre o instrumento e as paredes da caixa é de 
1,0 pol. Pergunta-se, se a caixa cai acidcníidrncnte de uma altura de 20 pol, o 
instrumento baterá nas paredes? 



Figura P.4-1-6. 


4-17 Um sistema mola-massa da Fig. P.4-17 tem um amortecimento de Coulomb que 
exerce um atrito constante dc força f. Mostrar que a solução para uma excita¬ 
ção dc ba.se c 

— = •—Lfl — COScu„í) — sena;„í ' 

('o cu„íiV mco/ 

onde a velocidade de base do Probl. 4-26 é admitida. 




Figura P.4-1 7. 
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4-18 Mostrar que a resposta pico para o Probl. 4-17 é 



Dividindo-se por a quantidade pode ser representada gra¬ 

ficamente como função de (^„ti com fti/mva como parâmetro. 

4-19 No Probl. 4-18, a força máxima transmitida a >n é 

-/-I- 


Representar graficamente esta quantidade cm forma não-dimensional, multipli¬ 
car por tijmvo para obter 


que novamonte pode ser representada graficamente como função de cor,, 
com parâmetros ftJmvQ. Traçar I cj„ ^niax/’’» I l^max/''oi'i 1 como fun¬ 
ção de aj„fi para/íi/mVfl igual a 0, 0,20 e 1,0. 

4-20 Mostrar que o espectro da resposta para o pulso retangular de duração í,, re¬ 
presentado na Fig. P.4-20 é dado por 


(•;í) 2scn^!". ív . 0,50 



4-21 O espectro da resposta para o pulso scnoidal é indicado na Fig. P.4-21. Mostrar 
que a resposta pico ocorre na região t > 1 1 para pequenos valores de í i /r. 
Mostrar que a resposta pico ocorre em t = t^, quando ti/r = 1/2. 



ü 12 3 4 5 6 


IJt 

Figura P.4-2J. 

4-22 Um sistema mola-massa não-amortccido com w' = 16,1 Ib tem um período 
natural de 0,5 segundos. Ele está sujeito a um impulso de 20 Ib/pol de forma 
triangular, com a duração de 0,40 segundo.s. Determinar o deslocamento máxi¬ 
mo da massa. 

4-23 No caso do um pulso triangular de duração ti, mostrar que, quando 



0 1 ' 2 3 . 4 5 6 


/r 

Figura P.4-23, 

t^lr = 1/2, a resposta pico ocorre em t = íi, que pode ser estabelecido a 
partir da equação 


1 






